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ABSTRACT 

We propose a new evolutionary model of a supermassive black hole (SMBH) and a circumnuclear 
disk (CND), taking into account the mass-supply from a host galaxy and the physical states of CND. 
In the model, two distinct accretion modes depending on gravitational stability of the CND play a key 
role on accreting gas to a SMBH. (i) If the CMD is gravitationally unstable, energy feedback from 
supernovae (SNe) supports a geometrically thick, turbulent gas disk. The accretion in this mode is 
dominated by turbulent viscosity, and it is significantly larger than that in the mode (ii), i.e., the CMD 
is supported by gas pressure. Once the gas supply from the host is stopped, the high accretion phase 
(~ 0.01 — O.lM0yr _1 ) changes to the low one (mode (ii), ~ lO _4 M0yr _1 ), but there is a delay with 
~ 10 8 yr. Through this evolution, the gas-rich CND turns into the gas poor stellar disk. We found 
that not all the gas supplied from the host galaxy to the central 100 pc region accrete onto the SMBH 
even in the high accretion phase (mode (i)), because the part of gas is used to form stars. Outflow from 
the circumnuclear region also suppresses the growth of the SMBH. As a result, the final SMBH mass 
(A/bh, final) is not proportional to the total gas mass supplied from the host galaxy (M sup ); MBH,final/Af sup 
decreases with M sup . This would indicate that it is difficult to form a SMBH with ~ 10 9 M© observed 
at high-z QSOs. The evolution of the SMBH and CND would be related to the evolutionary tracks of 
different type of AGNs. We found that the AGN luminosity tightly correlates with the luminosity of 
the nuclear starburst only in the high-accretion phase (mode (ii)). This implies that the AGN-star burst 
connection depends on the evolution of the AGN activity. 

Subject headings: galaxies:active — galaxies: nuclei — ISM:structure — galaxies:starburst — black hole 
physics 



1. INTRODUCTION 

The energy emitted by active galactic nuclei (AGNs) is 
commonly ascribed to accretion onto a supermassive black 
hole (SMBH). The recent discovery of powerful QSOs at 
z > 6 (Fan et al. 2001) implies that the formation of 
SMBHs is completed in less than 1 Gyr. In addition, 
high-resolution observations of galactic centers indicate 
the presence of a SMBH whose mass correlates with the 
mass of spheroidal stellar components of galaxies, i.e., 
Mbh/MspIi ~ 10~ 3 in the nearby universe (e.g., Kormendy 
& Richstone 1995; Marconi & Hunt 2003). These indi- 
cate that the formation and evolution of a SMBH and a 
spheroidal stellar component are closely related, despite 
nine orders of magnitude difference in their size scale. 

Various models have been proposed to explain the co- 
evolution of AGNs and spheroids, or the physical link be- 
tween the SMBH growth and the starburst in spheroids 
(e.g., Silk & Rees 1998; Adams et al. 2001; King 2003). 
However, little has been elucidated concerning the physics 
of angular momentum transfer in a spheroidal system (a 
bulge), which is inevitable for formation of SMBHs, tak- 
ing into account the tight connection between SMBHs and 
spheroids. Umemura (2001) proposed that the Mbh — 
M sp h relation can be explained by the mass accretion onto 
the galactic center from a galactic scale (~ 1 kpc) via 
the radiation drag (see also Kawakatu & Umemura 2002). 
During the hierarchical formation of a galactic bulge, the 



mass accretion due to the tidal torque driven by the ma- 
jor and minor merger of galaxies would also be an impor- 
tant process for SMBH formations (e.g., Toomre & Toomre 
1972; Mihos & Hernquist 1994; Mihos & Hernquist 1996; 
Saitoh & Wada 2004). However, previous studies of the co- 
evolution model simply postulated that the mass accretion 
onto the BH horizon, Mbh is regulated by the Eddington 
rate, MEdd (e.g., Umemura 2001; Kawakatu, Umemura, 
& Mori 2003; Granato et al. 2004), or equals the mass- 
supply rate from the host galaxy, M sup (e.g., Di Matteo et 
al. 2003; Kawata & Gibson 2005; Di Matteo et al. 2005; 
Springel et al. 2005; Okamoto et al. 2007). 

The gas accumulated by galaxy mergers and radiation 
drags does not accrete onto a SMBH directly, since the 
angular momentum of the gaseous matter cannot be thor- 
oughly removed. Thus, some residual angular momentum 
would terminate the radial infall, so the accreted gas forms 
a reservoir, i.e., a circumnuclear disk, in the central ^100 
pc around a SMBH (Fig. 1). In such a circumnuclear disk, 
the active star formation has been observed as the nuclear- 
starburst (< 100 pc) in nearby Seyfert galaxies (e.g., Iman- 
ishi & Wada 2004; Davies et al. 2007; Watabe, Kawakatu 
& Imanishi 2008). The nuclear starburst must affect the 
SMBH growth because radiation and/or supernova feed- 
back due to starbursts can trigger the mass accretion onto 
a SMBH (e.g., Umemura 1997; Wada & Norman 2002). 
In order to understand the SMBH growth in terms of the 
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evolutionary sequence of host galaxies, it is crucial to link 
mass accretion processes from a galactic scale with those 
from an accretion disk in the vicinity of a central BH, via 
the circumnuclear disk. 

To this aim, it is necessary to construct a model of the 
circumnuclear disk. Recently Thompson et al. 2005 (here- 
after T05) built up the radiation pressure-supported star- 
burst disk with AGN fueling, considering the mass supply 
from host galaxies. In their model, they postulate external 
torques (e.g., bars and spiral waves) in order to ensure that 
the gas can accrete to the galactic center before it turns 
into stars. We here propose another plausible model of 
a nuclear starburst disk supported by the turbulent pres- 
sure led by supernova explosions. Wada & Norman (2002; 
hereafter WN02) examined the structure of matter around 
AGNs associated with starburst by three-dimensional hy- 
drodynamic simulations. They showed that the internal 
gas motion is not steady and the matter is inhomoge- 
neous, but the global geometry is supported by internal 
turbulence caused by SN explosions. They also observed 
the mass accretion (~ O.lM yr _1 ) to the central a few 
pc region from the turbulent torus, and it is enhanced by 
the SN feedback. This means that the star formation and 
the gas accretion due to the turbulent viscosity can co- 
exist in this inhomogeneous disk. Our model presented 
here relies on this picture, in which angular momentum 
is transferred by a turbulent viscosity. In our model, the 
mass accretion process, star formation, and structure of 
the inhomogeneous circumnucler disk are naturally com- 
bined. We elucidate how SMBH grows from a seed BH, 
taking into account the mutual connection between the 
mass-supply from a host galaxy and the physical states 
of the circumnuclear disk accompanied by the star forma- 
tion. T05, on the other hand, did not discuss a long-term 
evolution of the circumnuclear disk and the central SMBH. 

This article is arranged as follows: In §2, we describe a 
new physical model for the circumnuclear disk supported 
by the turbulent pressure via SN explosions, considering 
the AGN fueling (§2.1, §2.2 and §2.3). In §2.4, we summa- 
rized the assumptions and free parameters in our model. 
In §3, we first show how the growth rate of a SMBH 
changes with time (§3.1). Then, we examine the evolution 
of SMBHs and the gas and stellar masses in the disk (§3.2). 
Next, we examine the relation between the final mass of 
SMBH and the total accreted gas mass from host galaxies 
(§3.3). In §3.4, we elucidate the time-evolution of AGN 
luminosity and the nuclear-starburst luminosity. By com- 
paring with observations, we predict the physical states 
of the circumnuclear disk for different types of AGNs. Fi- 
nally, we show the AGN-starburst luminosity relation, and 
then we discuss the origin of the AGN-nuclear starburst 
connection (§3.5). In §4, we discuss the effect of the AGN 
outflow and the radiation pressure caused by the young 
stars in disk. Moreover, we comment on the origin of a 
bias in AGN formation. Section 5 is devoted to our con- 
clusions. 

2. MODELS 

We presuppose that the dusty gas is supplied around a 
central SMBH at a rate of M sup from a host galaxy whose 
surface density, Shost, including the gas and stellar com- 
ponents, is constant in time. The accumulated gas forms 



a turbulent pressure-supported circumnuclear disk around 
a central SMBH with Mbh as schematically described in 
Fig. 1. In this paper, we focus on the average properties of 
the disk although we take into account the physics in the 
inhomogeneous distribution of gas in the disk as we will 
mention later. This treatment is valid because the global 
structure is dynamical stable even in the inhomogeneous 
distribution of gas in disk on a local scale (see WN02). 

2.1. Turbulent pressure-supported circumnuclear disk 

We describe the physical settings of the circumnuclear 
disk supported by the turbulent pressure via the SN ex- 
plosions. 

From the radial centrifugal balance, the angular velocity 
f2(r) is given by 



2 GM B h . 
= ; 1 (Sdisk(r) 



^host ) 



(1) 



where Mbh is the mass of the BH and HdiskM is the sur- 
face density of baryonic components (the gaseous matter 
and stars). 

Concerning the vertical structure of circumnuclear disk, 
we assume the hydrodynamical equilibrium. The turbu- 
lent pressure (Ptub = PgVt) associated with SN explosions 
is balanced to gravity, g caused by 



Ps( r ) v t( r ) = Ps(r)gh(r), 



(2) 



where v t (r) and h(r) are the turbulent velocity and the 
scale height of the disk, respectively. Here, the gravity, g 
is obtained as g = GM B nh/r 3 + 7rG(£ disk (r) + Ehost)- 

The geometrical thickness is determined by the balance 
between the turbulent energy dissipation and the energy 
input from SN explosions. Thus, the energy balance in 
unit time and volume can be obtained as 

Pg(f> t 2 (r) _ Pg(»> t 3 (r) 



= 7]S4r)E SN , 



(3) 



t d[s (r) h(r) 

where the dissipation timescale of the turbulence, ta s {r) 
is assumed to be a crossing time, h(r)/vt(r), £?sn is the 
total energy (10 51 erg) injected by an SN, n is heating ef- 
ficiency per unit mass which denotes how much energy 
from SNe is converted to kinetic energy of the matter, 
and the star formation rate per unit volume and time is 
S*(r) = C»p g (r). Here C* is the star formation efficiency 
and p g (r) = T, g (r)/2h(r) = f g T l< a s k(r)/2h(r) 7 where E g 
is the surface density of gas component in the disk and 
/ g is a gas fraction. In the inhomogeneous distribution 
of gas in disk, the star formation can be led by gravita- 
tional collapse of higher density clumps than the critical 
density. Recently, Wada & Norman (2007) showed that 
the star formation rate, C* is approximately proportional 
to the average density, p g especially for a high density end, 
if the statistical structure of gas density is represented by 
a log-normal probability distribution. 

By using eqs.(2) and (3), the turbulent velocity and the 
scale height for the region where central BH dominates the 
gravitational potential are given by 

w tM = 5— h tuh (r), (4) 



h t ub(r) 




(G^sn) 1/2 , 



(5) 
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where ht u b denotes the scale height supported by the tur- 
bulent pressure due to the SNe feedback. We should em- 
phasize that the validity of these solutions (eqs. (4) and 
(5)) were confirmed by the comparison with WN02. 
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Fig. 1. — Schematic picture for a model of SMBH growth linked 
with the change of physical state in the circumnuclear disk sup- 
ported by the turbulent pressure via SN explosions, where r; n and 
rout are the inner radius and the outer radius. Star symbols rep- 
resent the region where stars can form, i.e., r > r c , where r c is 
determined by Toomre criterion (see §2.2). Here, Mgjj, M g , Mbh, 
M, and M Bup are the mass of SMBHs, the gas mass, the stellar mass, 
the BH growth rate, the star formation rate, and the mass-supply 
rate from a host galaxy, respectively. 

2.2. Two regimes of gas accretion in circumnuclear disk 

We suppose a kinetic viscosity as a source of angular 
momentum transfer in the gas disk. Then, we adopt the 
formula of mass accretion rate in a viscous accretion disk 
(Pringle 1981) as 



M(r) = - 



lG{r,t) 



= 2TrvT, g (r) 



dlnO(r) 



d\nr 



(6) 



where the viscous torque generated by the turbulent mo- 
tion, G is defined as G(r, t) — 27w t (r)£ g (r)r 3 dil / 1 dr . Since 
the viscous parameter is expressed by v t (r) = avt(r)h(r), 
the mass accretion rate in general increases as the energy 
input from SN explosions becomes large, in other words, 
the star formation rate is high (see eqs. (4) and (5)). Thus, 
the mass accretion rate depends on the gravitational sta- 
bility of the disk. With respect to the stability, we adopt 
Toomre's stability criterion, i,e., when the surface density 
of gas in the disk, £ g is higher (lower) than the critical 
surface density, S cr ;t the disk is gravitationally unstable 
(stable). The critical surface density is obtained as 

E crit (r) = (7) 

7T(_r 

where n(r) = 4£l(r) 2 + 2f2(r)rdO(r) / dr is the epicyclic 
frequency and c s is the sound velocity. The critical ra- 
dius r c is determined by the Toomre criterion, that is, 
£ g (r c ) = S crit (r c ). 

In this picture, it is natural that there are two modes 
of gas accretion rate as follows: mode (i): If the circumnu- 
clear disk is fully gravitationally unstable, in other words, 



the critical radius for the gravitational instability, r c is 
smaller than the inner radius of the disk, n n , then the 
disk is geometrically thick due to stellar energy feedback, 
and as a result we have a large accretion rate. The value 
of a produced by randomly distributed SN explosions is 
not obvious. However, thanks to numerical simulations 
demonstrated by WN02, the average transfer of angular 
momentum can be described as v%(r) w v t (r)h(r). Thus, 
we suppose a = 0(1). mode (ii): If the critical radius is 

located inside the disk (r c > n n ), the scale height of the in- 
ner region would be much smaller than mode (i), because 
the scale height is determined by the thermal pressure, 
P g (r) = p g (r)ghth(f) where P g {r) — p g (r)c^. In the mode 
(ii), the magneto-rotational instability could be a source 
of turbulence, but the turbulent velocity is comparable or 
even smaller than the sound speed, i.e., a w 0.01 — 0.5 and 
v t = Cs (e.g., Balbus & Hawley 1991; Machida et al. 2000; 
Machida & Matumoto 2003). As a result, the accretion is 
less efficient than the mode (i). 

Since the critical radius and the inner edge of the disk 
are functions of the BH mass, and surface density of the 
gas, and the gas mass depends on the star formation rate 
in the gas disk, the evolution of the whole system (a central 
BH plus a circumnuclear disk) should be time-dependent. 

2.3. SMBH growth and states of the circumnuclear disk 

Based on the picture described in §2.1 and §2.2. we here 
examine the evolution of central BHs and circumnuclear 
disk, considering the mass supply from host galaxies (see 
Fig. 1). For this purpose, we evaluate the time-evolution 
of SMBH growth and the star formation rate and gas mass 
in the disk by focusing on the time dependence of charac- 
teristic radius in the disk, instead of solving the evolution 
of radial structure of disk. Then, the surface density of 
disk is assumed by a power law with the cylindrical radius 
r as Sdi 8 k(r) = S diski o(r/r out )"' ) ', where r out is the outer 
boundary of the disk, which will be defined below. Here- 
after the physical quantities with suffix are the value at 
Bout- 
in this model the supplied gas from the host galaxy is 
eventually consumed to form the SMBH or stars (see §4.1 
for the effect of AGN outflow). Then, the BH mass and 
gas mass in the disk are simply given by the mass conser- 
vation as follows: The time-evolution of the gas mass in 
the disk, M g = J r '° ut 27rr'£ g (/)dr', is given by 

M g (t) = / [M sup (i') - M*(t') - MbbWW, (8) 



where M sup (i), M*(i) and Mbh(£) are the mass-supply 
rate from hosts, the star formation rate, and the growth 
rate of SMBH, respectively. On the other hands, the time- 
evolution of SMBH mass Mbh (t) is obtained as 

M BH (t) =M BH ,scod+ / MBa{t')dt', (9) 
Jo 

where we assume the mass of seed BHs, Mbh, seed = 
10 3 M Q , as end-products of the first generation stars (e.g., 
Heger et al 2003). 

As for M sup , we can assume any function for the mass 
supply rate, but here we simply take a step function as the 
first attempt as 

const., for t < t s 



M sup (t) - i Q for ' t > t 



'-supj 



sup J 



(10) 
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where t sup is a period of the mass-supply from hosts. 

The star-formation rate in the disk can be defined as 
A/* = 2 J^° ut S*{r')h{r')2-Kr' dr' , where r min is the mini- 
mum radius at which stars can form. Here rmin = Hn for 
r c < r- m and r m ; n = r c for r c > r- m . Thus, M*ft) can be 
expressed in terms of M g (t) as follows: 

M,(i) = AC»M g (t), (11) 
where A — 1 for r c < r- m) while A = 1 — (r c /r out ) 3-7 for 
r c > r in . 

We define the growth rate of SMBH, M BH as M BH = 
M (rj n ), where Mftin) is a mass accretion rate in the inner 
region of the circumnuclear disk. By using eq. (6), the BH 
growth rate is calculated for < 7 < 2 as follows: 



gravitationally stable, and then no stars can form in the 
disk. Since our main aim is to reveal the physical con- 
nection between the AGN activity (the SMBH growth) 
and the starburst in the circumnuclear region, we con- 
centrate on inner parts of the circumnuclear disk (r < 
''out)- Then, when Mdisk > Mbh the outer radius r out is 
given by GM disk /r out = GM BB: /r out + TrGY, hos tr ou t, where 
J r out 27rr'E<3isk(''')'^ r ''j an( i r in is the inner radius 



disk 



(M disk / 7 rI] host ) 1 /2 



/ b 77£ S N(GM B H/rout)- 1 M*(n„/rout) 3 - 7 , for mode (i) 
I 3/ gco (Q^/G)(M g /M B H)(r in /r out ) 2 - 7 , 

where / b = 3(14 - 4 7 )(8 - 4 7 )/16 and / geo = 
2.4. 



M BH = 



0. 



for , mode (ii) 

(12) 
e r- m /h(r in ) 



of the disk. From eq. (2), r out is 

60pcAf^ )g E^^ 4 , where M dis k,8 and S hos t : 4 are the to- 
tal baryonic mass of the circumnuclear disk normalized 
by 10 8 M© and the surface density of host galaxy normal- 
ized by 10 4 Af Q pc~ 2 , respectively. On the other hand, 
if Mdisk < Mbh, ''out is obtained from GM B n/r out = 
7rG£host''out- From this, we quantitatively evaluate r out 
as r out = (Af B H/7rE hos t) 1/2 

Mbh, 8 is the mass of BH normalized by 10 s M©. 
Inner radius: r- m 

We suppose that the inner radius of the disk, n n , is deter- 
mined by the dust-sublimation radius, which is strongly 
supported by observational results (e.g., Suganuma et al. 
2006). Thus, the inner radius r- ln is expressed as ri n — 

T-L5QQ 8 , where Tagn,46 is the AGN luminosity 



60pcAf B g 8 £j/ t 2 4 , where 



1/2 



AGN, 46^1500 



Assumptions and free parameters 

In this section, we summarize the assumptions and free 
parameters as follows; 
Surface density of hosts: Ehost 

We fix the surface density of host galaxies with Shost = 
10 4 A/q pc~ 2 . This is the upper value of nearby starburst 
galaxies (e.g., Kennicutt et al. 1998), but may be a typi- 
cal value of high-z starburst galaxies (e.g., Tacconi et al. 
2006). 

Surface density of disk: Ediskft) 

As for the power law surface density profile of the disk, 
we should note that recent numerical simulations of AGN 
disks showed that the single power law profile is main- 
tained (e.g., Levine et al. 2007). Such a power low profile 
would be supported because the turbulence generated by 
the self-gravity and/or the SN feedback redistributes the 
angular momentum in the disk on a dynamical time scale. 
In this paper, we change the power law index 7 from to 
2. 

Gas temperature of disk: T g 

As for the gas temperature, we assume the isothermal gas 
with T g = 100 K, that is, c s cs lkms^ 1 . Recent nu- 
merical simulations showed that the circumnuclear disk 
around a SMBH would be dominated by the cold gas 
(T g ss 50 — 100 K) since the dust cooling is effective (e.g., 
Wada & Tomisaka 2005; see also Levine et al. 2007). In 
addition, the gas temperature is almost constant for the 
range of gas density in circumnuclear disks (Fig. 5 in Wada 
& Tomisaka 2005). These indicate that an isothermal gas 
with T f 

mass evolution of the circumnuclear disk 
Outer radius: r out 

The size r out is defined as the outer boundary inside which 
the potential of the BH plus circumnuclear disk system 
dominates that of the host galaxy. From Toomre con- 
dition, the outer part of the disk (r > r out ) is always 

2 Let us consider r; n , r ou t, Mg, M* and A/bh dependences on Mbh ( see e< 3-(12)). As smaller r ou t or larger r; n , S g (r; n ) increases because of 
S g (r) oc r -7 (0 < 7 < 2). The larger rj n (or smaller r out ) leads to larger Mjh, because Mbh ^ ^g{ r in) from eq. (6). The larger BH mass 
leads to a smaller scale height from eq. (5). Thus, Mbh becomes smaller since AfaH ex h. As M* (oc M g ) and M g are larger, E g (n n ) is larger, 
and thus A/bh increases. 



3 pcT 

in units of 10 46 ergs _1 and T1500 is the silicate grains sub- 
limation temperature in units of 1500 K (Laor & Draine 
1993). Although r m depends on the size and composition 
of dust grains (e.g., Barvainis 1992; Sitko et al. 1993; Laor 
& Draine 1993; Kishimoto et al. 2007), our conclusions 
does not change so much. Provided that £agn is equal to 
the Eddington luminosity, I/Edd = 47tcGA/bh't1p/0t, the 
inner radius can be rewritten as r; n = 3 pc Mg^gT^QQ 8 , 
where m p is the proton mass and ctt is the Thomson cross- 
section. In this paper, we adopt r; n = 3pcA/g H 2 g . 

In our model, we have three free parameters as follows: 
(i) the power law index of surface density, 7, (ii) the heat- 
ing efficiency, 77 and (iii) the star formation efficiency, C*. 
In the following sections, we use the fiducial value as 7 = 1, 
?y = IO^A/q 1 which can be derived from comparing our 
analytical solution of scale height with numerical simula- 
tions done by WN02 and C* = 3 x 10 -8 yr _1 which is upper 
value of nearby starburst galaxies (Kennicutt 1998; Fig. 13 
in Wada & Noman 2007). However, since 7, rj and C* are 
free parameters in this work we will discuss the effect of 
changing these values in §3.2. For these fiducial free pa- 
rameters, Mbh (eq. (12)) can be evaluated quantitatively 
as follows: 

0.3 77_ 3 C_ 8 E^ ti4 A/ g , 8 (t)A/ dis k,8(i)- 1/2 , 
for mode (i) 

-3„,„_„3 V V4 



g — 100 if is reasonable assumption for discussing 



Mbh ft) 
M Q yr" 1 



3 x 10- 3 a .5 <! E^ )4 M g)8 ft) A/ disk ,8ftr 1/2 
xMBH.sft) -1 ^ 4 , for mode (ii), 



where 77-3 is the heating efficiency normalized by 



10 3 M 1 , C* t -g is the star formation efficiency normal- 
ized by 10 _8 yr _1 , M gi8 is gas mass of massive disk nor- 
malized by 1O 8 M0, ao.5 is the a-parameter normalized by 
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0.5, and c S) i is the sound velocity normalized by 1 kms . 

3. RESULTS 

3.1. The growth rate of SMBHs and star formation rate 

Based on the physical model connecting between the 
growth of SMBHs and the physical states of circumnu- 
clear disks described in §2, we examine how the growth 
rate of SMBHs is related to the star formation rate in the 
circumnuclear disk. 

Figure 2 shows the time evolution of Mbh(£), M*(i), 
and M E dd{t) = L Edd (t)/c 2 for M sup = lMgyr 1 and 
isup = 10 8 yr. The critical time-scale, t a u is the time 
when r c — r m . Thus, the stars can form in the whole 
region of the disk before t = i C rit- We define the former 
phase (t < i C rit) as the high-accretion phase and the latter 
phase (t > t cr it) as the low-accretion phase. As seen in 
Fig. 2, we find that not all the gas supplied from a host 
galaxy reaches to the SMBH, i.e., M B h = O.l-O.3M yr _1 
for M sup = 1 AIq yr -1 even in the high accretion phase 
(r c < r; n ). This high accretion phase drastically changes 
to the low accretion phase with Mbh ~ lO _4 Af yr _1 . 
But, there is a delay with ~ 10 8 yr after the mass supply is 
stopped. The disk becomes gravitationally stable because 
the surface density of gas in the disk decreases due to the 
starformation. Thus, this timescale is basically determined 
by the star formation timescale, t* = C^ 1 ~ 3 x 10 7 yr. 
The drastic change of Mbh depends on whether stars are 
formed in the inner region of circumnuclear disk, because 
the ratio of Mbh between high and low accretion phase is 
given by 

M BH (for r c < r in ) V t h ub 



M BH (for r c > r in ) ac s h 



itr 



th 



where v t /(ac s ) « 10 2 M^ 8 and h tuh /h th 
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Fig. 2. — Time evolution of M B h (*), M*(t), M sup (t) and M Edd (t) 
lMgyr -1 and t Bup = 10 8 yr. t sup is a period of the 



stop 



for M s , 

mass-supply from hosts. i cr ; t is the time when r c 
the time when the massive disk is gravitationally stable, that is, 
S g < £ cr it for the entire massive disk. We call the phase (t < t c rit) 
the high-accretion phase, while we call the phase (t > t cr i t ) the low- 
accretion phase. The star formation efficiency C, is 3 X 10~ 8 yr. 

It also has been found that SMBH growth dominates 
the star formation in the disk (i.e., Mbh > M») in the 

Q " 1/2 

J During t < t Bup , Mbh M g because M^k is nearly equal to M g 
M g oc t until t = t BU p (see Fig. 3). 



early epoch, i.e., < < t» = C, 1 ~ 3 x 10 7 yr. In this phase, 
a super-Eddington mass accretion rate (Mbh 3> Medd) 
keeps for ~ 10 8 yr. Recent two-dimensional radiation hy- 
drodynamic simulations have shown that such a super- 
Eddington mass accretion flow is possible (e.g., Ohsuga et 
al. 2005). After this phase, M* is larger than Mbh because 
of M* oc t and M BH oc t 1 / 2 0. Then, the SMBH growth 
rate is much smaller than the Eddington mass accretion 
rate because of the gas pressure supported geometrically 
thin disk. After t = t s t op the star formation terminates be- 
cause the entire disk is gravitationally stable, where i s top 
is the time when S g (r) is smaller than E cr ; t (r) throughout 
the entire circumnuclear disk. 

3.2. Evolution of gas, SMBH and stellar masses 

In this section, we elucidate the physical difference be- 
tween two phases for SMBH growth (high- and low- ac- 
cretion phases). We plot the mass of SMBH, Mbh^), 
gas mass in the disk, M g (t), and stellar mass in the disk, 

M*(t) = f*M*(t')dt' in Fig. 3. The total supplied gas 
mass from the hosts, M sup = M sup t sup , is 1O 8 M . In 
the high accretion phase, plenty of gas accumulate around 
a SMBH since M sup (t) > Mbh(*) (see Fig. 2). Thus, 
the mass ratio of gas in circumnuclear disk and SMBH, 
mdisk = Mg/MaH is about ten and the gas fraction to the 
total baryonic mass, / g = M g /(M g + M„), is close to ~ 1. 
Recalling that the accretion is a super-Eddington in high 
accretion phase, large mdisk > 1 and / g ~ 1 would be 
the conditions which super-Eddington mass accretion can 
keep. Since the scale height in the inner region is inflated 
by the energy input led by the SN explosions, the cir- 
cumnuclear disk would be geometrically thick, which may 
correspond to the obscuring torus proposed by the AGN 
unified model (e.g., Antoucci 1993; Urry & Padovani 1995; 
WN02). On the other hand, in the low- accretion phase 
(t > t C iit), stellar mass dominates the BH mass and gas 
mass, i.e., / g « 10~ 2 and mdisk ~ 0.03. In this phase, the 
disk is the gravitationally stable, so that the inner region 
of disk should be geometrically thin. 

Moreover, we find that the saturation of SMBH growth 
appears ~ 10 s yr later after the mass supply is stopped. 
One should note that the final BH mass is about 1/3 of 
the total gas mass supplied from hosts. The timescale of 
the SMBH growth, £ grc .wth is determined by the viscous 



10 8 yr 



m 

v t h 



'BH, final 




10 8 M Q 

is the final SMBH mass. 

17/ 



VlO 8 yr, 
At given r\ 



lO^M^ 1 and M B H,finai = 3 x 10 7 M o (see Fig.3), the 



where M B H,finai 

timescale of SMBH growth is comparable to that of star 
formation in the gas disk. From eq. (14), it is found 
that the rapid BH growth (i.e., i grow th < 10 9 yr) can be 
achieved if the star formation rate in the disk is quite 
high, e.g., < 10 8 yr. Hence, it is essential to reveal what 
determines the star formation timescale (i.e., the star for- 
mation efficiency, C» ) in the circumnuclear disk. This will 

from eq. (12). On the other hand, M* oc M g from eq. (11), considering 
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be left as our future works. For the dependence of 77 and 
Mbh, final, we can understand as follows. From eq. (5), 
the larger BH mass leads to a smaller scale height, so that 
^growth becomes larger as Mbh, final increases because of 
Mbh oc h. Since Mbh oc 77 (see eq. (12)), £ gr0 wth is longer 
as 77 becomes larger. 

Finally, we discuss the effect of three free parameters (7, 
C* and rf). As for 7, we examine Mbh for 7 = 0.5 and 
7 = 1.5 by using the analytical solutions (see eq. (12)), 
and find that our results do not depend on the difference 
of 7 significantly, that is, the difference in Mbh between 
7=1 and 7 = 0.5 (or 1.5) is less than a factor of 2 although 
Mbh is larger as larger 7. For C*, both the star formation 
rate and and BH growth rate in the high accretion phase 
are ten times smaller from eqs. (11) and (12), if we adopt 
ten time smaller C», i.e., C, = 3 x 10~ 9 yr _1 . Thus, after 
t sup the gas mass in the disk decreases slower (see eq. (8)), 
and then /j cr ;t and t s top become roughly ten times longer. 
However, the final BH mass (Mbh, final ~ MsHigrowth) 
do not change significantly because of Mbh ck C* and 
^growth oc C~ . This is confirmed by Fig. 4, which is the 
same as Fig. 2 and Fig. 3 but the star formation efficiency 
is ten times smaller (i.e., C, = 3 x 10~ 9 yr _1 ). Concerning 
77, if 77 is ten times smaller than the value we adopted, i.e., 
10" 4 Mq 1 , the final BH mass is ten times smaller because 

of Mbh oc V in the high accretion phase (see eq. (12)). 
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Fig. 3. — Time evolution of the mass of BH (black solid line), 
-MbhMi the gas mass in the disk (blue solid line), M g (t), and the 
stellar mass in the disk (red solid line), M*(t). The total sup- 
plied mass from host galaxies M aup is 1O 8 M0. The star forma- 
tion efficiency C* is 3 X 10~ 8 yr. Note that m^isk = Mg/M^H and 
f g = M e /(M g + M,). 
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Time [yr] 

Fig. 4. — Same as Fig. 2 and Fig. 3, but the star formation 
efficiency C* is ten times smaller (i.e., C, = 3x 10 — 9 yr _1 ). 

3.3. M B H,finai-io-M sup relation 

As mentioned in §3.2, the final mass of SMBH, MBH,finab 
is not equal to M sup , i.e., Mbh, final ~ 0.3M sup . In this 
section, we examine the relation between Mbh, final and 
M sup . Figure 5 shows the final BH mass as a function 
of total accreted mass from host galaxies with M sup = 
10 6 M Q -10 10 M Q . We here assume t sup = 10 8 yr. Since the 
solid line is our prediction and the dot-dashed line denotes 
M B h, final = M sup , it is found that M B H,finai/M sup gets 
smaller as M sup becomes larger, i.e., Mbh final /M sup = 
0.8,0.6,0.3,0.1 and 0.05 for M sup = 10 6 ,'l0 7 , 10 8 , 10 9 , 
and 10 10 M Q , respectively. In other words, the final mass of 
SMBHs, Mbh, final is not proportional to M sup and growth 
of a SMBH is more inefficient for larger M sup . As men- 
tioned in §3.1 and §3.2, the star formation in the disk 
dominates the growth of the SMBH as M sup increases. 
As a result, all the gas accreted from hosts cannot ac- 
crete onto a central SMBH, but turns into stars in the late 
phase of mass-supply. The above trend does not change 
for t sup = 10 7 — 10 9 yr and also for changing the star for- 
mation efficiency, C*. 

The averaged mass-supply rate, M sup , is expected to 
be typically w 0.1 — 1-Mgyr , driven by galaxy merg- 
ers (e.g., Mihos & Hernquist 1996; Saitoh & Wada 2004) 
or by the radiation drag associated with the starbursts in 
hosts (e.g., Umemura 2001; Kawakatu & Umemura 2002). 
Thus, it may be hard to accrete the gas mass more than 
w 10 9 M Q from the host galaxy for 10 8 — 10 9 yr. Even if the 
gas supply with lM Q yr _1 lasts for 10 9 yr, a SMBH does 
not evolve more than 10 8 M Q . This would imply that it is 
difficult for the only gas accretion process due to the tur- 
bulent viscosity driven by SN explosions to form a SMBH 
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> 10 9 M Q observed in the luminous high-z QSOs (e.g., 
McLure et al. 2006). In other words, the direct gas ac- 
cretion from a host galaxy into the accretion disk (sub-pc 
scale) and/or the mergers of compact objects for the for- 
mation of SMBH with > 1O 8 M . 
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Fig. 5.— The final SMBH mass, Mbh, final against the dif- 
ferent total accreted mass from host galaxies, Af sup = M sup t sup 
for t Bup = 10 8 yr. The thick solid line shows our prediction. The 
dot-dashed line denotes Mbh final = M SU p. We show the parameter 
space of M sup and Mbh final which are hard to achieve by any plausi- 
ble mass-supply process (M sup > 1O 9 M0 and Mbh, final > 1O 8 M0), 
as the shaded region. Note that our prediction does not depend on 
C» and tsup (see §3.2 and 3.3). 



3.4. Evolution of AGN luminosity and nuclear starburst 

luminosity 

We examine the evolution of the AGN luminosity and 
the nuclear-starburst luminosity (=the luminosity origi- 
nating in star formation in the circumnuclear disk), in 
order to understand which objects correspond to these 
two phases observationally. In this paper, we define these 
luminosities as follows: The AGN luminosity (=the ac- 
cretion disk luminosity) can be obtained as a function 
of tobh = Mbh/Mem (Watarai et al. 2000), by as- 
suming two types of accretion disk, namely the slim disk 
(Abramowicz et al. 1988) and the standard disk (Shakura 
& Sunyaev 1973). 



' 2(l + ln™^W)i Edd (t) ;m BH (t)>20, 



£agn(*) = < 

L (*^)Wt) ;™ BH (t)<20. 

(15) 

The nuclear starburst luminosity Lnuc Sb(£) can be given 
by 

iNuc,SB(t) - 0.14eM»(t)c 2 , (16) 

where e = 0.007 which is the energy conversion efficiency 
of nuclear fusion from hydrogen to helium. 

Figure 6 shows the evolution of AGN luminosity 
£agn(£) and nuclear-starburst £nuc,sb(£) for M sup = 

The mass accretion in units 



yr. 



lMgyr 1 and t sup = 10 s 
of Eddington mass accretion rate itibh is larger than 20 
until t = t sup (see Fig. 2). Then, Lagn(£) depends 
on the evolution of LEdd(i) °c Mbh(£) from eq. (15), 
and thus Lagn(^) increases with time. During this phase 



(t < t S up), the nuclear-starburst luminosity, Lnuc,Sb(£) oc 
M*(t) oc M g (t), also increases with time because M g (t) 
increases (see Fig. 3). After t — t sup , both Lagn(£) and 
£nuc,sb(£) decrease since M g starts to decrease monoton- 
ically (see Fig. 3). The AGN luminosity decreases drasti- 
cally after t = t CI i t , reflecting the phase change of accre- 
tion (see Fig. 2). After t = t stop £nuc,sb becomes zero 
because the circumnuclear disk is gravitationally stable 
as mentioned in §3.1. In summary, there are four phases 
in Lagn(£) and -£nuc,sb(£); (i) the super-Eddington lu- 
minosity phase: Lagn > ^Edd {t < t sup ), (ii) the sub- 
Eddington luminosity phases: Lagn = (0.1 — l)^Edd 
(tsup < t < t C rit), (ni) the starburst-luminosity dominated 
phase: Lagn < Ln uc ,sb < ^Edd (*crit < t < t stop ), and 
(iv) the post AGN/starburst phase: Lagn *C ^Edd & 

Lnuc,SB = (t > tstop)- 
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Fig. 6. — Time evolution of the AGN luminosity, LagnM, that 
of the nuclear starburst luminosity, Lnuc Sb(*) an d that of the Ed- 
dington luminosity, ^EddW- We can divide this into the four phases 
as follows: (i) t < t sup , (ii) tsup < t < t clit , (iii) t cr it < t < tstop 
and (iv) t > t s top- These phases (i), (ii), (iii) and (iv) correspond 
to the super-Eddington luminosity phase, the sub-Eddington lumi- 
nosity phase, the starburst-luminosity dominated phases and the 
post- AGN/Star burst phase, respectively. 

Phase (i) might be a good model for narrow-line type 
I Seyfert galaxies (NLSls) because NLSls are super- 
Eddington AGNs (e.g., Mathur ct al. 2001; Kawaguchi 
2003). On the other hands, if the mass supply is extremely 
high (e.g., M SU p = 10M Q yr _1 ) and consequently the cir- 
cumnuclcar region is covered by dense gas, then this phase 
might be observed as ULIRGs. Note that some kind of 
ULIRGs show properties similar to NLSls (e.g., Teng et 
al. 2005; Hao et al. 2005; Kawakatu et al. 2007a). If 
this is the case, we predict that a gas-rich massive torus 
(wdisk > 1 and / g = 1) exists at their center. Such a mas- 
sive torus can be observed through carbon monoxide and 
hydrogen cyanide molecular emission using the Atacama 
Large Millimeter Array instrument (ALMA) for NLSls 
and ULIRGs (e.g., Wada & Tomisaka 2005; Kawakatu et 
al. 2007b; Yamada, Wada & Tomisaka 2007). Judging 
from the luminosity ratio (iAGN/^Edd = 0.1 — 1), in phase 
(ii) the AGN is still bright and thus this might be BLSls 
or QSOs. In this phase, the relatively small mdisk and / g 
are predicted, i.e., mdisk = 0.2 — 1 and / g w 0.1 — 0.5. 
After phase (ii), the AGN luminosity becomes smaller by 
three order of magnitude than its peak luminosity. Thus, 
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phases (iii) and (iv) may correspond to the low luminosity 
AGNs (LLAGNs) associated with and without the nuclear 
starburst, respectively. As seen in Fig. 3, a geometrically 
thin massive stellar disk, i.e., M* > Mbh and f s « 10 -2 
is expected in this post-AGN phase. Thus, this indicates 
that the LLAGNs do not have obscuring tori envisaged in 
the AGN unified model because no stars form near the in- 
ner radius of the circumnuclear disk. In other words, we 
predict the absence of obscuring tori in LLAGNs. This 
seems to be consistent with observations that most low lu- 
minosity radio galaxies (FR I) show no evidence of dusty 
tori expected in normal AGNs (e.g., Chiaberge et al.1999). 
As an alternative explanation, Elizur & Shlosman (2006) 
suggested the torus disappears when the AGN luminosity 
decreases below ~ 10 42 ergs _1 because the mass accre- 
tion in the accretion disk cannot maintain the required 
cloud outflow (see also Honig & Beckert 2007). In or- 
der to understand the origin of absence of dusty tori in 
LLAGNs, it is essential to explore the spatial distribution 
of young stars in them. In addition, since such a stellar 
disk is a remnant of a gas-rich massive torus in the early 
phase of SMBH growth, its discovery would strongly sup- 
port the idea that LLAGNs are dead Seyfert galaxies or 
QSOs. Therefore, the observations of the physical states 
of a circumnuclear disk on < 100 pc, e.g., Wdisk and / g , 
for different type of AGNs would be useful to understand 
AGN evolution. 

Finally, we summarize a coevolution scenario of SMBHs 
and circumnuclear disks in Fig. 7. According to our 
knowledge obtained in this paper, we might guess that 
the final state of AGNs depends on the mass supply rate 
M sup and the duration of the mass-supply t sup . In details, 
we will discuss this issue in our on-coming paper. 




Fig. 7. — Schematic sketch for the coevolution of a SMBH and 
a circumnuclear disk. Star symbols (yellow) represent young stars. 
Star symbols (red) denote old stars. The darkness of gray color 
shows the gas fraction (/ g ), that is, the darker gray represents higher 
gas fraction. 

3.5. AGN luminosity vs. Nuclear- starburst luminosity 

relation 

Based on our model, we present the correlation between 
AGN luminosity and nuclear-starburst luminosity. In Fig. 
8, we plot £nuc,SB against Lagn for the different Mbh, final 
(or M sup ). The higher concentration of data points rep- 
resents that the evolution is slow. We find that AGN lu- 
minosity positively correlates with nuclear-starburst lumi- 



nosity for bright AGNs (right branches of Fig. 8), which 
corresponds to the high-accretion phase, that is, phase (i) 
and phase (ii) defined in §3.4. This trend is consistent 
with the observational results for nearby Seyfert galaxies 
(Imanishi & Wada 2004; Watabe et al. 2008). Moreover, 
-^NucjSb/ Lagn becomes larger as Lagn increases. This 
can be explained as follows: Lnuc,sb/Lagn (oc M*/Mbh) 

1 /2 

is larger as M g increases because of M*/Mbh oc M g ' (see 

1 /2 

§3.1). Considering Lagn oc Mbh oc M g ' during t < i sup , 
-^Nuc,sb/£agn increases with Lagn- We should mention 
that it is worth examing this trend by future observations. 

On the other hand, the Lagn vs - ^Nuc,sb relation in 
LLAGNs (left branches) does not follow that in bright 
AGNs. In this phase, the range of AGN luminosity is 
independent of Ln uc ,sb, i-e., Lagn = 10 41 — 10 42 ergs _1 
for Lnuc.sb = 10 40 — 10 44 ergs _1 , because A/bh is basi- 
cally determined by the sound velocity in the disk. Since 
the nuclear-starburst luminosity, Z.x„ c .sb oc M* oc M g in- 
creases with M g , no tight correlation between Ln U c,sb and 
Lagn is seen in LLAGNs. In summary, we suggest that 
the AGN-starburst connection depends on the evolution of 
the AGN activity. In order to test our model, it is essen- 
tial to compare the AGN-starburst luminosity relation for 
bright AGNs (e.g., Seyfert galaxies and QSOs) with that 
for LLAGNs. 
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Fig. 8. — The nuclear starburst luminosity, Lnuc,SB against the 
AGN luminosity, Lagn- The dots with different colors denote the 
different final BH mass, A/bh, final- The dashed lines represent the 
luminosity ratio, Lnuc Sb/^AGN = 0.01, 0.1 and 1, respectively. 

4. DISCUSSION 

4.1. Effects of AGN outflow 

Through this paper, we assume all the gas is accreted by 
a central BH, i.e., Mbh = M(n n ) even when the mass ac- 
cretion rate at r m exceeds the Eddington limit (We call this 
model "model (a)"). According to the recent two dimen- 
sional radiation hydrodynamic simulations (e.g., Ohsuga 
et al. 2005), the super-Eddington mass accretion actually 
can be achieved, but they mentioned that the large part 
of accreting gas turns into outflow because of the strong 
radiation pressure from the accretion disk (Ohsuga et al. 
2005; Ohsuga 2007). 

We here assume that the dynamics of circumnuclear disk 
is not affected by the outflow. The two models of AGN 
outflow are possible when M(n n ) > MEdd as follows; 
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(i) The BH grows by a super Eddington rate but the 
part of accreting matter forms outflow, "mode/ (bf . 

Outflow = fy,M(r in ) and 

M BH = (1 - U)M(r in ), 
where we assume f w = 0.1 from the results of Ohsuga 
(2007). 

(ii) The growth of BHs is regulated by the Eddington 
limit and the remnant of accreting gas turns into the out- 
flow, "model (cf. 

Moutflow = M(r in ) - Mbh and 

Mbh = CgjljMEdd, 

where cbh = 0.05 is the energy conversion efficiency. 

Figure 9 is the same as Fig. 2 and Fig. 3, but the out- 
flow rate (M ou tfl ow ) and the total outflow mass (M out fl w) 
are considered based on model (b), where M out fl ow = 
J * MoutRow{t')dt' . During t < t of the SMBH grows by 
the super-Eddington mass accretion rate, but most of ac- 
creting mass from the gas disk is blown away as the AGN 
outflows, where i f is the time when M(r- m ) = c^MEdd- 
After t = t f , all the gas accrete onto a central BH by less 
than the Eddington limit (i.e., no AGN outflow). After 
t = tcrit the evolutions are the same as the case without the 
outflow, i.e., model (a). As seen in Fig. 9, we find that the 
final mass of SMBHs is Mbh, final — 8 x 10 6 M Q , which is a 
factor four smaller compared with the prediction by model 




(a). 



Fig. 9. — Same as Fig. 2 and Fig. 3, but the outflow rate 
(Moutflow) an d the total outflow mass (M out fl ow ) are newly added 
for the model (b), i.e., super-Eddington growth of SMBHs with the 
AGN outflow. Here i f is the time when M(n n ) = e^A^Edd- 

In the model (a), the super-Eddington accretion mainly 
determines the final SMBH mass (see Fig. 2 and Fig. 3). 



But, in model (b) the contribution of super-Eddington ac- 
cretion is comparable to that of Eddington accretion for 
the final SMBH mass. This indicates that the contribu- 
tion of super-Eddington accretion depends on the strength 
of AGN outflows, i.e., / w . Therefore, it will be worth re- 
vealing which physics controls / w by using the radiation 
hydrodynamic simulations. 

On the other hands, we show the case of Eddington lim- 
ited BH growth (i.e., model (cf) in Fig. 10. Apparently, 
the BH growth is greatly suppressed by AGN outflow and 
consequently the final SMBH is order of magnitude less 
than model (a), namely A/bh. final — 1O 6 M0. In summary, 
we find that the evolution of SMBHs are closely linked 
with not only the state of circumnuclear disk, but also the 
strength of AGN outflow and the conditions under which 
the super-Eddington growth occurs. 




Fig. 10. — Same as Fig. 2 and Fig. 3, but the outflow rate 
(Af ou tfl ow ) and the total outflow mass (M ou tfl w) are newly added 
for the model (c), i.e., Eddington limited growth of SMBHs with 
the AGN outflow. 

4.2. Effects on the UV radiation pressure 

So far, we have ignored the effect of radiation pressure, 
but the UV radiation from young stars in the gas disks 
may affect on the vertical structure of disk. Thus, we ex- 
amine whether the radiation pressure is really negligible 
when the disk geometry is determined by the SN feed- 
back. If the disk is optically thick to the UV radiation, 
the radiation pressure (P ra d) is given by 

Pad = Isb/c ~ 0.14e5,(r)ftt ub (r)c, (17) 
where h tu b(r) is determined by eq. (5). Here Isb = 
0.14eS*(r)h(r)c 2 is the starburst luminosity per unit area. 

As a result, the ratio of the radiation pressure and the 
turbulent pressure via SN explosions (Ptub = Pg v t) can De 
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obtained as follows: 

^^:W^y /2 , (X8) 

where idyn = {r 3 /GMbk) 1 ^ 2 is the dynamical timescale. 
Because of idyn(?"*out)A* ~ 10~ 2 for the relatively high star 
formation efficiency, i.e., C* = 10 _8 yr , the radiation 
pressure is neglected as far as r\ > IO^Mq 1 . This condi- 
tion (rj > IO^Mq 1 ) is satisfied for a turbulent pressure 
supported starburst disk modeled by WN02 using three- 
dimensional hydrodynamical simulations. Thus, our con- 
clusions do not change significantly even if we consider the 
UV radiation pressure. 

4.3. Origin of a bias in AGN formation 

Thanks to a growing number of observations at differ- 
ent redshifts in multi- wavelength, it has been suggested 
that the density of more luminous AGNs associated to 
more massive SMBHs peaks earlier in the history of the 
universe than that of low luminosity ones (Miyaji, et al. 
2000; Ueda et al. 2003; Hasinger et al. 2005; Bongiorno et 
al. 2007) . This evolution has been dubbed as AGN cosmic 
downsizing. It is still unclear why the growth timescale of 
SMBHs is shorter in more massive SMBHs, in other words, 
why the growth time is highly biased in an environment 
to form massive SMBHs. 

From eq. (14), the final SMBH mass can be described 

as 




The AGN downsizing might be explained by our 
result, if larger Mbh, final grows in the circumnuclcar 
disk with a higher star formation efficiency. For in- 
stance, Fig. 11 shows that the growth timescale for 
-A/bh. final = lO 8 A/0 is shorter than that for Mbh, final = 
1O 7 M , i.e., ig r owh(M B H,finai = 1O 8 M ) m 10 8 yr and 

Wh(^BH,final = 10 7 M Q ) M 10 9 yr, if C,(M B H,final = 

1O 8 M ) is w 30 times as large as C*(M BH ,finai = 10 7 M Q ). 
According to Wada & Norman (2007), they showed that 
the star formation rate increases extensively as a function 
of average gas density. Since the massive galaxies could 
be formed in an environment with a higher average gas 
density, the star formation efficiency could be 10 times 
larger if the average gas density is ss 10 times higher. As a 
consequence, the UV radiation field would be stronger in 
more massive galaxies. The recent radiation hydrodynam- 
ics simulations showed that the critical density for the star 
formation can be affected by the UV radiation (e.g., Susa 
& Umemura 2004). These works indicate that the critical 
density is larger for a stronger UV radiation field. In addi- 
tion, the star formation efficiency, C* would be larger for a 
larger critical density at given average gas density (Wada 
& Norman 2007). Therefore, since the circumnuclcar disk 
in more massive galaxies is exposed to the stronger UV 
radiation, the star formation efficiency is larger for the 
gas disk in massive galaxies. Therefore, the strength of 
UV radiation field surrounding the galactic nuclei might 
control the timescale for SMBH growth. In order to clar- 
ify this issue, it is essential to examine the star formation 
process in the gas disk under an intense UV field with the 
three-dimensional radiation transfer. 
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Fig. 11. — The final SMBH mass, Mbh final against the timescale 
of the SMBH growth without the AGN outflow (model (a)). The 
solid lines with different colors correspond to the different star for- 
mation efficiency. The shaded region denotes the parameter space 
of AfflH final which is difficult to achieve in the present model. 

5. CONCLUSIONS 

We have elucidated the physical role of the circumnu- 
clear disk in relation to SMBH growth. To this end, we 
have constructed a new model of SMBH growth, taking 
into account the evolution of physical states of the cir- 
cumnuclear disk formed by the mass-supply from the host 
galaxy. In the present model, we consider the two regimes 
of gas accretion depending on the gravitational stability of 
the disk. Our main conclusions are summarized as follows: 

1. Not all the gas in the disk accrete onto the SMBH, 
i.e., M B n = 0.1-0.3M Q yr" 1 for M sup = 1 M yr" 1 
because the part of gas is used to form stars in 
the disk. This high accretion phase changes to the 
low accretion phase with Mbh ~ 10 -4 M Q yr. But 
the transition takes ~ 10 8 yr after the mass sup- 
ply from hosts is stopped. This timescale is basi- 
cally determined by the star formation timescale, 
t* ~ C^ 1 ~3x 10 7 yr. Through this evolution, the 
gas rich disk turns into the gas poor stellar disk. 
In the high accretion phase, a super Eddington ac- 
cretion is possible, and thus the existence of gas 
rich circumnuclcar disks (M g > Mbh) is a condi- 
tion which the super Eddington accretion onto a 
central SMBH keeps. In the present model, the two 
phases of SMBH growth depend on whether stars 
can form in the inner region of the circumnuclear 
disk. This would imply that the AGN activity de- 
pends on the spatial distribution of young stars in 
the circumnuclear disk. In order to examine this, in- 
frared observations with high spatial resolution (e.g. 
Very Large Telescope Interferometer (VLTI) , Space 
Infrared Telescope for Cosmology & Astrophysics 
(SPICA) and 3m Atlantic infrared telescope) are 
crucial, because young stars would be buried in the 
optically thick circumnuclcar disk. 

2. Evolution of nuclear activities are divided into four 
phases; (i) the super-Eddington luminosity phase, 
(ii) the sub-Eddington luminosity phase, (iii) the 
starburst-luminosity dominated phases, and (iv) 
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the post-AGN phase. Judging from the luminos- 
ity ratio, Lagn / LEdd, NLSls or ULIRGs may cor- 
respond to phase (i). If this is the case, it is pre- 
dicted that they have gas-rich tori, whose masses 
are larger than the mass of their central SMBHs. 
It indicates that the SMBH growth proceeds in the 
phase surrounded by gas-rich tori. Phase (ii) may 
correspond to BLSls or QSOs. In this phase, we 
predict that the gas fraction of the circumnuclear 
disk is smaller than that for NLSls or ULIRGs, i.e., 
/ g ~ 10 — 50%. Finally, phases (hi) and (iv) may 
correspond to the low-luminosity AGNs (LLAGNs). 
According to our prediction, LLAGNs have gas- 
poor (gas fraction f g ~ 1%) circumnuclear disks. 
This result suggests that the physical states of cir- 
cumnuclear gas or stellar systems are important to 
understand the nature of different types of AGNs 
for gas components with ALMA and/or VLTI for 
stellar components. 

3. The hnal SMBH mass (Mbh, final) is not propor- 
tional to the total gas mass supplied from the 
host galaxy (M sup ) during the hierarchical forma- 
tion of galaxies, i.e., MBH,finai/Mmp decreases with 
M sup , because the star- formation rate overcomes 
the growth rate of SMBHs as M sup increases. Con- 
sidering the plausible mass accretion processes from 
the host galaxy, our results may indicate that it is 
difficult to form w 10 9 M Q SMBHs observed at high- 



z QSOs, only by the gas accretion process due to the 
turbulent viscosity via SN explosions. This might 
imply the importance of alternative mechanisms for 
the formation of - 10 9 M Q SMBH. 

4. The AGN luminosity (Lagn) tightly correlates with 
the nuclear-starburst luminosity (£nuc,sb) in the 
bright AGNs (Seyfert galaxies and QSOs). We also 
predict that £n uc ,sb/£agn is larger for more lumi- 
nous AGNs (i.e., QSOs). This is not the case in the 
LLAGNs. Therefore, the preset model predicts the 
£agn vs. Lnuc,sb relation depends on the evolution 
of AGN activity. 

5. The AGN outflow from the circumnuclear disk 
greatly suppresses the growth of the SMBHs, es- 
pecially if the BH growth is limited by the Edding- 
ton rate when the mass accretion from the disk ex- 
ceeds the Eddington limit. This indicates that the 
strength of AGN outflows affects on the evolution 
of SMBHs and hnal mass of SMBHs. 
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